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Abstract 

We study propagation dynamics of a particle phase in a single-file pore connected to a 
reservoir of particles (bulk liquid phase). We show that the total mass M(t) of particles 
entering the pore up to time t grows as M(t) = 2m(J, pp)y/Dot, where D is the "bare" 
diffusion coefficient and the prefactor m(J, pp) is a non-trivial function of the reservoir 
density pp and the amplitude J of attractive particle-particle interactions. Behavior of the 
dynamic density profiles is also discussed. 

Introduction 

Particles transport across microscopic pores is an important step in a vast variety of 
biological, chemical engineering and industrial processes, including drug release, catalyst 
preparation and operation, separation technologies, especially biological and biochemical, 
tertiary oil recovery, drying and chromatography |], [| [| . 

Man-made or naturally occuring porous materials contain a wide range of pore sizes, 
from meso- to micro- or even nanoscales, in which case the pore diameter is comparable to 
the molecular size. Such molecular sized channels, of order of a few Angstroms only, appear, 
for instance, in biological membranes, and are specific to water and ion transport which 
participate in hydrostatic or osmotic pressure controlled cellular volume regulation [|J. 
Carbon nanotubes or zeolites, such as, Mordenite, L, AIPO4-5, ZSM-12, may also contain 
many channels of nearly molecular diameter and can selectively absorb fluids serving as 
remarkable molecular sieves |J. For example, AIPO4-5 is composed of nonintersecting and 
approximately cylindrical pores of nominal diameter 7.3 Angstroms. 

A salient feature of transport in molecularly sized pores is that there is a dramatic 
difference between the diffusion of adsorbates whose size is much smaller than the pore 
diameter, and those whose size is comparable to it. If the diameter of the diffusing guest 
molecules exceeds the pore radii, but is still less than its diameter, the particles are able 
to enter the pore but not able to bypass each other such that initial given order is striktly 
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Figure 1: Solid with a single- file pore in contact with a liquid phase. 



maintained (see Fig.l). Here, diffusion obeys the so-called " single- file" behavior; the mean- 
square displacements of tracer molecules do follow X 2 (t) = 2Ft 1 ^ 2 || 0, where F is 
referred to as the "single- file" mobility; evidently, such a behavior is remarkably different 
from the conventional diffusive law X 2 (t) = 2Dt, observed in situations when the guest 
molecules are able to bypass each other. Single-file diffusion was evidenced by pulsed 
field gradient NMR measurements with a variety of zeolites and guest species || [7j; in 
particular, the law X 2 {t) = 2Ft 1 / 2 has been identified experimentally for C2H 6 || and CF 4 
@ in AIPO4-5. 

Diffusion of absorbed particles in single-file pores has also been the issue of a consider- 
able theoretical interest recently. Several approaches have been proposed, based mostly on 
the lattice-gas-type models, and such properties as concentration profiles or steady-state 
particle currents have been evaluated [TDl O, O. As well, a great deal of Monte Carlo and 



Molecular Dynamics simulations has been devoted to the problem |13|, [TJ|, |Tj| , providing a 
deeper understanding of the transport mechanisms in single-file pores. On the other hand, 
still little is known about non-stationary behavior in single-file systems; in particular, how 
fast do the particle phase propagates within the single-file pores or how does the total mass 
(or number) of particles entering the pore up to time t grows with time? 

In the present paper we focus on the challenging question of the particle phase propaga- 
tion dynamics in the single-file pores. We report here some preliminary results; a detailed 
account will be published elsewhere ||16|| . More specifically, we consider a single-file pore in 
contact with a reservoir of particles (bulk liquid phase), which maintains a fixed particle 
density at the entrance to the pore. The pore is modelled, in a usual fashion, as a one- 



dimensional regular lattice whose sites support, at most, a single occupancy; the particles 
interaction potential consists of an abrupt, hard-core repulsive part, which insures single- 
occupancy, and is attractive, with an amplitude J > 0, for the nearest-neighboring particles 
only. Introducing then a standard, interacting lattice-gas dynamic rules (see e.g. fT7|), we 
derive evolution equation for the local variables describing mean occupation (density) of 
the lattice sites, which is analysed both analytically and numerically. Within our appoach, 
we define the evolution of the total number M(t) of particles entering the single-file pore 
up to time t and also discuss the dynamics of the density distribution function within the 
pore. We show that the growth of M(t) is described by M(t) = 2m( J, pp)^D t, which 
can be thought of as the microscopic analog of the Washburn equation, where D Q denotes 
the "bare" diffusion coefficient, while the pref actor m(J, pp) is a non-trivial function of the 
attractive interactions amplitude J and the density pp at the entrance to the pore. 



The model 

Following earlier works 0, [LI], [12| , as well as a conceptually close analysis of an upward 
creep dynamics of ultrathin liquid films in the capillary rise geometries JL8|| , which appears 



to be well-adapted to the single-file dynamics, we model the single-file system under study 
(Fig.l) as a semi-infinite linear chain of equidistantly placed sites X (with spacing a), 
attached to a reservoir of particles maintained at a constant chemical potential p (see 
Fig. 2). Each pair of sites is separated by a potential barrier of height Eg, which sets the 
typical time scale tb- Note that the spacing a can be defined, in case of hard solids, as the 
interwell distance of a periodic potential describing the interactions of the particles with 
the solid atoms and tb is related to E B and the reciprocal temperature (3 = 1/T through 
the Arrhenius formula. For soft solids, in which case the dominant dissipation channel is 
due to mutual particle-particle interactions, a can be thought of as the typical distance 
travelled by particles before successive collisions; here, t b is just the ballistic travel time. 

Reservoir 
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Figure 2: Effective model for transport in a single-file pore. 
Further on, we define the particles interaction potential U(X) as: 



0, X > a, 

U[X): { —J, X = a, (1) 
■foo, X < a, 

i.e. we suppose that the interaction potential between particle is a hard-core exclusion, 
which prevents multiple occupancy of any site, and attraction with an amplitude J, (J > 
0), between the nearest-neighboring particles only. Occupation of the site X at time 
moment t for a given realization of the process will be described then by the Boolean 
variable T] t (X), such that 

, f 1, the site X is occupied, /nS 

^ (X) = \0, otherwise. (2) 

Consequently, the interaction energy U t (X) of the particle occupying at time t the site 
X for a given realization of the dynamical process is 

U t (X) = -JU{X + a)+ Vt (X- a)). (3) 



Lastly, we define the particle dynamics (see |17j for more details). We suppose that at 
time moment t any particle occupying site X waits an exponential time with mean tb and 
then selects a jump direction with the probabilty 

p(X\X') = Z- 1 e W [^(u t (X)-U t (X'))}, $>(X|X') = 1, (4) 
i x , 

where Z is the normalization, X'{= X±a) denotes here the target neighboring site and the 
sum over X' means the sum over all nearest neighbors of the site X. As soon as the target 
site is chosen, the particle attempts to hop onto it; the hop is instantaneously fulfilled 
if the target site is empty; otherwise, the particle remains at its position. Physically, it 
means that repulsive interactions are very short-ranged - much shorter than the lattice 
spacing, and particles "learn" about them only when they attempt to land onto some 
already occupied site. In turn, attractive interactions are felt within the distance equal 
to the lattice spacing and hence, influence the choice of the jump direction in order to 
minimize the total energy of the system. 



Evolution equations 



Now, let pt{X) = r] t (X), where the overbar denotes averaging with respect to different 
realizations of the process. Assuming local equilibrium, we find then that the time evolution 
of Pt(X) is governed by the following balance equation fly| : 



T B f>t{X) 



1 - p t {X)) p t {X - a)p{X - a\a) + p t (X + a)p(X + a\X) 



Pt (X) (1 - Pt (X + a))p(X\X + a) + (1 - p t (X - a))p(X|X - a) , (5) 



where the average transition rate p(X\X') obeys Eq.([|) with r) t (X) replaced by pt(X). 
Equation @ accounts for the fact that a particle may appear at time moment t on an 
empty site X by hopping from the occupied sites X ± cr with corresponding transition 
probabilities dependent on the interaction energy of the system; and may leave the occupied 
site X for unoccupied sites X ± a. 

Next, we turn to the so-called diffusion limit, assuming that tb scales as a 2 ; we thus 
suppose that o — > 0, tb — > 0, but the ratio ct 2 /tb = const = 2Dq, where Z?o is the diffusion 
coefficient describing motion of an individual, isolated particle. Expanding pt(X ± a) and 
p(X|X ± a) in the Taylor series up to the second order in powers of the lattice spacing a, 
we arrive at the desired dynamical equation of the form 



p t (x) = A 



_d_ 

dX 



A- + Pptim - Pt (x))^l 



where U t (X) is the interaction energy at point X defined by Eq.fl3|). 

Now, a few comments on Eq.@ are in order. Note first that Eq.( 
equation with an environment dependent force, 



(6) 



is a Burgers-type 



dU t (X) 



-13 



dpt(X) 



-a 2 3 



d 2 U t (X) 



dX " dX " " dX 2 

When only the first term on the rhs of Eq.(^) is taken into acount, we get from Eq. 
a one- dimensional diffusion equation 



(7) 



Pt(X) 

with a field-dependent diffusion coefficient 



^-D(pt(X))^-pt(X), 
dX VPV JJ dX HK h 



D( Pt (X)) = D (l - 2p3 Pt (X)(l-p t (X))), (9) 

which is precisely the equation derived earlier by Lebowitz et al [19| and describing hy- 
drodynamic limit dynamics of a system of mutual interacting particles undergoing ballistic 
motion. On the other hand, if we keep the second term on the rhs of Eq.(0), (which is 
appropriate if we consider some steady-state solutions fIE[|), we will obtain the customary 
equation of the form 



(10) 



where the mobility M(p t (X)) is given by 



M(fH(X))= Pt (X)(l-p t {X)), (11) 
while local free energy J-(p) obeys 

a 2 (3Jsdp^ 



Hp) = Jdx(f( P ) + a -fL(§) 2 ), (12) 



with 



f(p) = p hip + (1 - p) In (1 - p) + (3Jp{\ - p). (13) 

Curiously enough, f(p), which has been derived in our work starting from a microscopic 
dynamical model obeying the detailed balance condition, has exactly the same form as the 



phenomenological Flory-Huggins-de Gennes local free energy density |2T[. Note also 
that for (3 J > 2, the local free energy f(p) in Eq. (|T3"D has a double-well structure whose 
minima approach and 1 as (3 J increases. This implies that the Onsager mobility in 
Eq.(|TTD never reaches negative values, contrary to the behavior predicted by Eq.(^) for 
which one has D(p) < when p c _ < p < p Cj+ , with 

Pc± = \(l±y/l-2/PJ). (14) 

Finally, we define the appropriate boundary conditions. As a matter of fact, any parti- 
cle, in order to enter to the nanopore from the liquid phase, has to surmount an additional 
barrier Em related to the enthalpic energy difference between the particle within the pore 



and in the bulk liquid phase [11, Supposing that the reservoir (bulk liquid phase) is in 



equilibrium with the particle phase within the single-file pore, we thus stipulate, following 



a similar analysis in ||18|| , that the reservoir maintains a constant density pp (see |L1], [12] for 
relation between pp and energetic parameters) at the entrance of the pore (site X = in 
Fig. 2). Second boundary condition is rather evident, we just suppose that Pt(X) vanishes 
as X — > oo at fixed t. Consequently, Eq.(|B|) (or Eq.flTOD) is to be solved subject to the 
conditions 

p t (X = 0)=p F , Pt (X^oo) = 0. (15) 
Below we discuss solutions of Eqs.(^) and ( |TUD obeying these two boundary conditions. 

Results 



We focus here on the time-evolution of the total mass of particles M(t), having entered 
the single-file pore up to time t, i.e., 



M(t) = / dX p t {X) (16) 

JO 

To define the time-dependence of M(t), it is expedient to turn to the scaled variable 
to = X/2y/D()t. In terms of this variable Eq.fllOD attains the form 



d 2 p{uj) dpiyj) d r / \dp{uj 



du 



^4D t' do; 



I du; 3 



while the boundary conditions in Eq.fllSD become p{uj = 0) = Pf and p(u — > oo) = 0. 

Note now that the term in the second line on the rhs of Eg. QT7D , associated with the 
second term in the expansion of the interaction energy in the Taylor series, Eq.([7|), is 
irrelevant to the dynamics as t — > oo, since it is multiplied by a vanishing function of 
time. Hence, the propagation dynamics can be adequately described by Eq.flD, which 
is not sufficient, however, for description of the steady-state characteristics, such as, e.g. 
steady-state particle current through a finite pore [[16] . 



Now, in terms of the scaled variable u, the total mass of particles M(t) reads 



M(t) = 2m(J,p F )y/D t, (18) 
where the prefactor m(J, pp) is determined by 

/•CO 

m(J,pF)= I duj p(u) (19) 
Jo 

Consequently, Eq.fllBD, which can be thought of as the microscopic analog of the Wash- 
burn equation, signifies that the mass of particles grows in proportion to the square-root 



of time. Note that similar result has been obtained in [IS| for J = 0, in which case 
m(J,p F ) = pp/Vn. 

Before we discuss behavior of the prefactor m(J, pp), it might be instructive to under- 
stand what are the physical processes underlying the M(t) ~ \/t behavior (see also 



for more detailed discussion). To do this, let us first recollect that the boundary condition 
Pt(X = 0) = pp = const is tantamount to the assumption that the reservoir is in equilib- 
rium with the particle phase in the pore. Turning next to the model depicted in Fig. 2, we 
notice that a jump of the rightmost particle of the particle phase away of the reservoir, 
leads to creation of a "vacancy". When this vacancy manages to reach diffusively, due to 
redistribution of particles, the entrance of the pore, it perturbs the equilibrium and gets 
filled by a particle from the reservoir. Hence, the mass of particles within the pore M(t) 



is proportional to the current J of vacancies from the front of the propagating phaseQ, 
M(t) ~ J7", where the current J ~ 1/L, L being the distance travelled by the rightmost 
particle away from the entrance of the pore. Consequently, M(t) ~ L ~ 1/L, which yields 
eventually the M(t) ~ \fi law. Note also that, from the viewpoint of the underlying 
physics, the dynamical process under study is intrinsically related to such phenomena as 
directional solidification, freezing, limited by diffusive motion of the latent heat, or Stefan 
problem. 

To define the prefactor m(J, pp) we have solved Eq. (|l7D (with the last term set equal 
to zero) numerically, for various values of the system's parameters pp and J. These results 
are summarized in Fig. 3 where we depict m(J,pp) as a function of pp for several different 
values of (3 J. 

Now, one notices that for any values of pp and (3 J the prefactor m(J, p F ) > 0, which 
implies that, as one can expect on intuitive grounds, there is no transition in the one- 
dimensional system under study and the particle phase propagates into the pore as soon 
as pp > 0. On the other hand, the prefactor m(J, pp) depends on the system's parameters 
in a quite non-trivial fashion. While for relatively small f3 J the prefactor m(J,pp) varies 
with pp almost linearly, for large f3J some saturation effect occurs, followed by, for larger 
(3 J, a non-monotoneous /^-dependence, and eventually, after a cusp-like variation, with a 
rapid growth. 

We begin with a small (3 J limit, in which case some perturbative analytical calculations 
are possible. To do this, let us represent p{uj) in the form of the series 

oo 

pH = ^(2/3J)> n H, (20) 

n=0 

(which is related to expansion in powers of the reservoir density pp < 1, as we will see in 
what follows) and try to calculate explicitly several first terms in such an expansion, con- 
straining ourselves to the quadratic in f3J approximation. After some rather cumbersome 
but straightforward calculations, we find eventually, that in the quadratic with respect to 
the parameter 2(3 J approximation, the prefactor m(J, pp) obeys: 



(2 p Jf 



(21) 



m(J,p F ) = - {2(3 J) [0.18/4 - 0.134 pp 

0.025 p 3 F - 0.047 p 4 p + 0.018 p% 

1 Note, however, that there are some subtleties concerning propagation of the rightmost particle of the 
phase growing in the single-file pore. As a matter of fact, it has been shown in |l8| that in a similar 
model without attractive interactions (i.e. J = 0) the mean displacement of the rightmost particle follows 
X(t) ~ yjt In (t), i.e. grows at a faster rate than the mass M(t). Similar behavior is observed for the front 
of the pahse, propagating in the single- file pore, also in the interacting case, i.e. for arbitrary J fig]. 




Figure 3: The prefactor m = m(J, pp) versus the density pp for different values of the amplitude 
J of the attractive particle-particle interactions. Symbols denote the results of numerical solution 
of Eq . (fl7|) , while the solid lines correspond to the analytical result in Eq . (|2l|) . The inset displays 
the behavior for small pp. 



This dependence, as one notices, agrees quite well with the numerical solution for 
relatively low values of (3 J over the entire domain of variation of pp, or, for larger (3 J, for 
progressively smaller values of pp. 

To get some understanding of the intricate, non-trivial behavior of m(J,pp), observed 
for larger values of (3 J, we analyse dynamical density profiles defined by Eq.® versus X 
and X/M(t) for (3 J = 3 and different reservoir densities pp, (see Fig. 4). Note now that 
the form of the density profiles is rather complex and depends largely on whether pp is 
less than or exceeds p c ,±. When pp < p c ,-, the form of the density profile is well described 
by Pi(X) = pi?erfc(X/M(t)), where erfc(X/M(t)) is the error function. This behavior is 



essentially the same as the one predicted for non-interacting lattice gas in fll8[ . On the other 
hand, when p F exceeds p C; _, but is less than p c |+ (for (3 J = 3, p Cj ± are equal to 0.79 and 0.21, 
respectively) we have two different regimes: the density rapidly, within the small constant 
distance 1(pf), drops to the value p Cj _ and then evolves as Pt(X) = p c _erfc(X/M(t)), where 



p Cj _ is independent of pf- Since, the prefactor m(J, pf) is just an integral of pt(X), we 
have, hence, two contributions: the first is the integral over the interval [0,l(pp)}, which 
is weakly dependent on pp and the second one - the integral over [l(pp), oo[, which is 
independent of pp and gives the bulk contribution to m(J, pp). These two contributions 
define the plateau-like part in the dependence of m(J, pp) on pp. On the other hand, when 
Pf exceeds p Cj+ and hence D(p) is positive definite, a dense droplet-like structure emerges 
near the entrance of the pore, which grows in size in proportion to \fi and contains most 
of the particles. This phenomenon explains apparently the growth of m(J, pp) with pp 
observed for p F > p c + . 

Conclusions 

In conclusion, we have studied propagation dynamics of a particle phase emerging in a 
single-file pore connected to a reservoir of particles (bulk liquid phase). Modelling the pore 
as a semi-infinite one-dimensional regular lattice, whose sites support, at most, a single 
occupancy, and supposing that the particles interaction potential consists of an abrupt, 
hard-core repulsive part, and is attractive, with an amplitude J > 0, for the nearest- 
neighboring particles, we derived the dynamical equation describing the time-evolution 
for the local particle density. This equation has been analysed both analytically and 
numerically. Within our appoach, we defined the evolution of the total number M(t) 
of particles entering the single-file pore up to time t and also discussed the dynamics of 
the density distribution function within the pore. We have shown that the growth of 
M(t) is described by M(t) = 2m( J, pp)y/D t, which can be thought of as the microscopic 
analog of the Washburn equation, where D denotes the "bare" diffusion coefficient, while 
the prefactor m(J,pp) is a non-trivial function of the attractive interactions amplitude J 
and the density pp at the entrance to the pore. We have discussed a peculiar behavior 
of m(J, p F ) and explained it through the analysis of the dynamical density distribution 
within the pore. 
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Figure 4: Dynamical density profiles defined by Eq.(||) for (3 J = 3 and different values of pp. 
In the left column pt{X) is plotted versus the space variable X and here different curves on each 
graph show the dynamics of the density distribution in the single-file pore, Eq.(||), while in the 
right column we present corresponding plots for pt(X) vs the scaling variable X/M(t). 



